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ABSTRACT 

S: 

, The Space Interferometry Mission (SIM) will make precise astrometric measurements that 

\ can be used to detect planets around nearby stars. We have simulated SIM observations and 

estimated the ability of SIM to detect planets with given masses and orbital periods and mea- 
' sure their orbital elements. We combine these findings with an estimate of the mass and period 

distribution of planets determined from radial velocity surveys to predict the number and char- 
■ acteristics of planets SIM would likely find. Our predictions are based on extrapolating the 

mass distribution of known extrasolar planets by up to a factor of ~ 100. This extrapolation 
provides the best prediction we can make of the actual number of planets that SIM will detect 
and characterize, but may substantially over- or underestimate the frequency of Earth-mass 
planets, especially if these form by a different mechanism than giant planets. We find that SIM 
j—{ | is likely to detect ~ 1 — 5 planets with masses < 3M ffi (depending on mission parameters). SIM 

would measure masses and orbits with 30% accuracy for ~ — 2 of these planets, but is unlikely 
to measure orbits with 10% accuracy for more than one of them. SIM is likely to detect ^5 — 25 
planets with mass less than 20Af®, measure masses and orbits with 30% accuracy for ~ 2 — 12 
of these, and measure masses and orbits with 10% accuracy for ~ 2 — 8 such planets. SIM is 
likely to find ~ 25 — 160 planets of all masses, depending on the observing strategy and mission 
lifetime. Roughly 25 — 65% of the planets detected by SIM have sufficiently large masses and 
short orbital periods that they can also be detected by radial velocity surveys. Radial velocity 
surveys could measure orbital parameters (not including inclination) for 30 — 70% of the planets 
whose orbital parameters will be determined to within 30% by SIM. 

Subject headings: planetary systems - astrometry - techniques: radial velocities 
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1. Introduction 

High precision radial velocity surveys have discovered ~ 100 planets around nearby main- 
sequence stars (Butler et al. 2002 and references therein; http://exoplanets.org), providing numer- 
ous challenges to theories of planet formation and evolution. The radial velocity technique is best 
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suited for massive planets in short-period orbits around lower main-sequence stars. Astrometric 
observations offer an alternative means for detecting extrasolar planets around nearby stars. Since 
the amplitude of the astrometric perturbation by a planet increases with the planet's orbital ra- 
dius and period, astrometric searches are more sensitive to long-period planets than radial velocity 
surveys if their time baseline is long enough. 

The Space Interferometry Mission (SIM) is expected to make high precision targeted observa- 
tions that will permit the detection of low mass planets around nearby stars. Indeed, the detection 
of low mass planets is the focus of two of the SIM key projects (Marcy et al. 2002; Shao et al. 2002) 
and one of the main science drivers for SIM. Since SIM observing time will be limited and very 
valuable, it is important to examine possible observing programs and their likely outcomes. In this 
paper, we estimate the number and properties of planets that SIM is likely to detect for several 
choices of mission parameters. 

While SIM is expected to achieve 4 /was absolute astrometry with wide-angle observations, it is 
expected that narrow-angle observations with ~ 1 fias differential astrometry will be most valuable 
for planet searches. In the narrow-angle mode, SIM observations will be grouped into roughly one 
hour intervals during which the baseline of the SIM interferometers will be maintained by locking 
onto bright guide stars. During each one hour interval it will be possible to make up to ~ 10 
one dimensional measurements of the angular separations of a target star and four reference stars. 
Photon noise is smaller than systematic errors for stars with V < 10.5. For pairs of stars with 
V < 10.5, the precision of an individual 30 second observation is expected to be < 3 fias. We 
define the single measurement precision to be the precision achieved by combining multiple nearly 
simultaneous observations in the usual way, assuming that the observational errors are random and 
uncorrelated. Thus, if ten measurements in a one-hour interval are devoted to a single target star, 
then the single- measurement precision is expected to be ~ 3fias/y/T0 ~ 1/uas. Given the numerous 
demands for SIM observing time, it is expected that the SIM planet searches will collectively be able 
to search ~120, 240, or 480 stars with a 1, 1.4, or 2 fxas single measurement precision, respectively 
(http:/ / sim.jpl.nasa.gov/ao_support/index.html) . 

Sozzetti et al. (2002) have simulated SIM observations to determine the probability that SIM 
would detect a planet given its orbital period and the amplitude of the astrometric perturbation it 
induces on its parent star. We use their results to determine the ability of SIM to detect planets 
around nearby stars, given the actual distances and masses of stars that SIM is likely to survey. 
We randomly assign planets to these stars based on the observed frequency of extrasolar planets. 
The masses and orbital periods of the simulated planets are drawn from the planetary mass and 
period function determined by Tabachnik & Tremaine (2002) based on radial velocity observations 
of massive planets. 

We describe our assumptions and simulations in §2. In §3 we present the results of our 
simulations. Based on these results we answer several questions regarding the discovery potential 
of SIM in §4. In §5 we summarize our main findings and conclusions. 
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2. Methods 

We have conducted Monte Carlo simulations to anticipate the number and properties of ex- 
trasolar planets that the narrow angle SIM planet search is likely to detect. First, we generate a 
list of target stars appropriate for SIM. Next, we assign hypothetical planets, specified by a mass 
and orbital period or semi-major axis, to these stars using an empirical fit to the known extrasolar 
planet mass-period distribution from Tabachnik & Tremaine (2002). Then, we determine which of 
these planets SIM would discover using empirical fits to the detailed simulations of Sozzetti et al. 
(2002). Finally, we analyze the number and the properties of planets discovered to anticipate the 
discovery potential of SIM. For the sake of comparison, we apply a simplistic model to estimate 
the number and properties of planets which could be discovered by a radial velocity survey of the 
same stars. 

2.1. Target Stars 

We begin by selecting stars from the Hipparcos catalog within 100 pc of the Sun (ESA 1997). 
Next, we restrict the sample to stars brighter than V = 10.5. For dimmer stars the photon noise 
would dominate the instrumental uncertainties (~ 1.7 //as) in a 30 second integration, demanding 
longer integration times to achieve the specified astrometric precision. Provided that V < 10.5, 
there is little gain in astrometric precision or observing efficiency by observing brighter stars, due 
to systematic uncertainties (http://sim.jpl.nasa.gov/ao_support/index.html). We require that the 
stars not be a member of a known spectroscopic binary or have an observed companion within 1". 
Of the stars which meet the above criteria, we select iV/4 stars of each of the four spectral types F, 
G, K, and M, for a total of N stars in the sample. Since the Hipparcos catalog includes luminosity 
classifications for only about half of these stars, we do not exclude giant stars. We estimate that 
~ 5%, and not more than 10% of the stars in our samples are giant stars. The actual SIM target 
list would likely replace these stars with main-sequence stars at slightly larger distances. 

Since obtaining a single measurement precision of ad < 3 //as requires averaging the results of 
multiple 30 second observing blocks, the necessary observing time scales approximately as ~ a^ 2 
for ad < 3 /(/as. Assuming constant mission time, the number of stars that can be surveyed therefore 
scales as ~ a d , provided that there is no additional overhead (e.g., slewing, acquiring guide stars). 
While a target star and its reference stars must fit within a ~ 1° diameter field of view, several 
target stars can be observed across a ~ 15° diameter field of regard (FoR) without slewing or 
significant additional mission time. Within a one hour observing interval, the angles between one 
target star and four reference stars (all V < 10.5) may be measured at a precision of cr^ = 1 //as. 
Alternatively, in the same hour interval two or four target stars could be measured at a precision 
of ad = 1.4 or 2 //as, respectively, provided that all the target and reference stars lie within a single 
FoR (http://sim.jpl.nasa.gov/ao_support/index.html). 

Thus, we take N =120, 240, and 480 for ad =1, 1.4, and 2 //as, respectively. For the latter 
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two samples, we must also ensure that one or three additional target stars are within each FoR. 
For these target lists we first choose closest candidate stars (the "primary stars") and then choose 
the nearest two or four stars that satisfy the other criteria above and fit in the required FoRs (the 
"secondary stars"). This requirement accounts for the tails of distant stars for the 1.4 ^uas and 2 
//as target lists shown in Fig. 1. The Hipparcos catalog is incomplete at magnitudes > 9, but this 
incompleteness is only beginning to become significant for the M stars in the larger target lists. 
We will consider a five-year mission with 24 two dimensional observations per target star and a 
ten-year mission with 48 observations per target star. 

We also consider an observing strategy similar to that of the Extrasolar Planet Interferometric 
Survey (EPIcS) proposed by Shao et al. (2002). This proposal uses a two-tier approach to observe 
~ 75 of the nearest stars (tier I) at 1 /uas precision 35 times, as well as ~ 2100 nearby stars (tier II) 
at 4 /was precision 25 times (over a five-year mission). The tier II stars all share an FoR with one 
of the tier I stars, making more efficient use of SIM observation time. We consider a similar, but 
simplified two-tier target list. To construct a survey that is directly comparable to the single-tier 
surveys discussed above, we choose 52 of the nearest stars (tier I) to observe at 1 uas precision and 
1092 nearby stars (tier II) at 4 ^xas. Each set of 21 tier II stars must be located within the 15° 
diameter FoR centered on a tier I star in our simulations. We observe both the tier I and tier II stars 
24 times during a five-year mission or 48 times during a ten-year mission. The one-tier and two-tier 
surveys which we consider require the same total number of individual 30 second integrations (i.e., 
120 ~ 52 + 1092/16). While the EPIcS survey plans to choose tier II targets to ensure a significant 
number of several types of stars, for our simulations of a two-tier survey we do not impose such 
constraints. In particular, our tier II target list includes 49 A stars, 213 F stars, 377 G stars, 354 
K stars, and 99 M stars. 



2.2. Planet Model 

Radial velocity surveys have discovered ~ 100 extrasolar planets (Butler et al. 2002) and 
measured their semi-major axis, a, and the ratio of the mass of the planet to the mass of the 
parent star times the sine of the orbital inclination of the planet relative to the sky, rasini/M*, 
where M* is the stellar mass. The values of msini for these planets range from 0.12Mj ~ 40M e 
to over 10Mj, where Mq, Mj, and Mq are the masses of the Sun, Jupiter, and the Earth. There 
appears to be a sharp decrease in the frequency of planets more massive than ~ 10Mj. Tabachnik 
& Tremaine (2002) have fit the planet mass ratio-period function with a power law distribution, 

dn = 0.0018 ( m \ ~°' 12 ( £-) °' 26 ^ (1) 
V0.0018M*y \8fld/ m P p v ' 

with dn = for orbital periods P p < 2 d or mass ratios m/M± > 10~ 2 ~ 10Mj/M Q . We 
extrapolate this planet mass function to the smaller masses that will be detectable by SIM. We 
use this functional form to assign planets, their mass ratios, and their orbital periods to the target 



- 5 - 



stars. A planet's semi-major axis is computed from its orbital period using Kepler's third law, 
assuming a stellar mass based on the spectral type (Hansen & Kawaler 1994). 

Note that Tabachnik & Tremaine (2002) assumed that the stars surveyed for planets all had 
mass equal to one solar mass. Thus, an alternative distribution that is equally consistent with their 
analysis is 

where M* in equation (1) has been replaced by M Q . If this alternative planet mass function were 
correct, then the number of planets of all masses that SIM would detect and characterize would 
be slightly increased, but the number of detectable planets with masses less than 20M e would be 
significantly decreased. 



2.3. SIM Model 

We estimate the ability of SIM to detect a planetary companion following the simulations of 
Sozzetti et al. (2002). They define a, the semi-amplitude of the star's astrometric wobble, by 

" — AU D ' U 

where m is the mass of the planet, M± is the mass of the parent star, a p is the semi-major axis of 
the planet, and D is the distance to the system from the observer. 

Sozzetti et al. (2002) present their results in terms of the "scaled signal" , 

S=-, (4) 

where o& is the single measurement precision of a set of one dimensional relative delay measurements 
between the target star and each reference star. Sozzetti et al. (2002) assume ~ 2 jtxas, but their 
results can be scaled for other values of ad- Assuming 24 two-dimensional observations equally 
spaced over the mission and an orbital period shorter than the mission duration, they found that 
a 50% probability of detecting a planet requires S > 1 and a 95% probability of detecting a planet 
requires S > 2.2. The Sozzetti et al. (2002) simulations detect a planet when the best fit model 
with no planets can be rejected by a x 2 test with 95% confidence. Since the target lists that we 
consider include 120 to 1144 stars, a 5% false alarm rate is unacceptably high. We prefer to raise 
the detection threshold of Sozzetti et al. (2002) by a factor 1.5, to S > 3.3, which gives a negligible 
false alarm rate (~ 0.1%), according to our own simulations. 

According to the simulations of Sozzetti et al. (2002), estimating the mass and orbital param- 
eters to within 30% of the actual parameters for 95% of the planets requires S > 5. To measure 
the mass of a planet accurate to within 10% with 95% probability requires S > 22. To measure 
the mass and orbital parameters to within 10% with 95% probability requires S > 33. 
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We have conducted our own simulations, which verify the main findings of Sozzetti et al. 
(2002). We agree that S > 33 is required to measure the mass and orbital parameters to within 
10% with 95% probability. We also find that the mass and orbital parameters can be measured 
to within 10% with 50% probability when 5^6. Thus, for many planets there is a significant 
chance that their orbits will be measured, even for S significantly smaller than the thresholds given 
by Sozzetti et al. (2002). We use our own simulations to estimate the probability of detecting or 
measuring the mass and/or orbital parameters for a given signal S (Ford 2003). We approximate 
the probability of detecting a planet or measuring its mass or orbital parameters based on a fit 
with the functional form 

P=\ I 

Ll + (£?/&)-/» + (S/c) -7 

where b, c, (3, and 7 are parameters obtained from fitting to the results of our own simulations 
and are listed in table 1. Different values of these fit parameters are used depending on whether 
the probability is for detecting a planet, measuring its mass and orbital parameters with < 30% 
accuracy, measuring its mass with < 10% accuracy, or measuring its mass and orbital parameters 
with < 10% accuracy. 

For orbital periods comparable to or longer than the mission lifetime (Pml), the astrometric 
signature required for a detection or measurement increases. To account for this effect we have 
divided our simulations according to the ratio of the planet's orbital period, P p , to Pml (see Table 
1). We use different values of b, c, /3, and 7 for each category. While it is clearly possible to 
detect the astrometric wobble induced by a planet with an orbital period much longer than the 
observational baseline, our simulations show that it is difficult to measure the orbital parameters 
if the observations do not span at least a quarter of an orbital period. Thus, we have imposed a 
sharp upper limit that prevents measuring orbital parameters — but allows detections — of any 
planets with orbital periods more than four times the mission lifetime. We allow detections for 
orbital periods up to twelve times the orbital period; the exact value of this cutoff has little effect 
on the results. 

We use the above formalism to calculate the probability that a planet will be detected or 
have its mass and orbital parameters measured, given the planet mass (m) and orbital period (P p ). 
Below we will present results for several sets of mission parameters. We consider single measurement 
precisions of aa = 1, 1.4, and 2 fias for target lists of 120, 240, and 480 stars, respectively. We 
consider a standard mission lifetime Pml = 5 yr with 24 two dimensional observations of each star 
as well as an extended mission lifetime Pml = 10 yr with 48 two dimensional observations of each 
star. We also consider the two-tier strategy described at the end of §2.1. For the ten- year mission 
we reduce the values of ad by a factor of \/2- Our simulations show that this closely approximates 
the effect of doubling the number of observations from 24 to 48. 



(5) 
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2.4. Radial Velocity Model 

To estimate the extent to which planets detected by SIM will be distinct from those found by 
radial velocity searches, we simulate the results of a radial velocity survey for each of our SIM target 
lists. Although precision radial velocity surveys of thousands of nearby stars are already underway 
at several ground-based observatories, often with 50 or more observations per star, for simplicity we 
assume the same number of observations for radial velocity surveys as for SIM, and that the radial 
velocity observations span the same time interval as SIM measurements. Due to these assumptions, 
we expect that we underestimate the number of planets which would be detected by actual radial 
velocity surveys. This is in contrast to the SIM surveys, for which we have chosen our assumptions 
to err on the side of overestimating the number of planets detected. 

To estimate which of the planets detected by SIM would likely be detected by radial velocity 
surveys, we apply a model similar to that used for the SIM astrometric survey, but with the scaled 
signature being given by 

where K is the star's velocity semi-amplitude, and <jk — 3m s" 1 is the accuracy of a single radial 
velocity measurement. We use the same functions to estimate the probability of detecting a planet 
or measuring its mass and orbital parameters as we use for SIM. (Strictly, radial velocities arc 
sensitive to the combinations msini and uj + rather than the mass, m, inclination, i, argument 
of pericenter, u, and longitude of ascending node, £1, individually.) We find the use of equation (6) 
gives reasonable agreement with the simulations of Cummings et al. (2002) assuming 24 observations 
over five years and a 0.1% false alarm rate, and are also roughly consistent with the current 
radial velocity discoveries (Butler et al. 2002 and references therein). The factor of \/2 in the 
denominator of the scaled signal for radial velocity observations accounts for the fact that radial 
velocity observations are only sensitive to one component of the star's velocity, while astrometric 
observations measure the star's two dimensional position. For the ten-year mission we assume 48 
observations and reduce the value of <tk by \/2. 

While this analysis neglects many important factors that affect the efficiency of radial velocity 
surveys, it allows us to estimate crudely what fraction of planets found by SIM could be indepen- 
dently discovered by a radial velocity survey of the same stars. Of course, SIM will allow the mass 
and inclination to be measured separately which cannot be done with radial velocity measurements 
alone. We do not consider the benefits of combining radial velocity and astrometric observations 
(Eisner & Kulkarni 2002). 

2.5. Summary 

Throughout this paper we attempt to make generous assumptions that tend to overestimate the 
number of detections and measurements SIM will make. For example, the simulations of Sozzetti 
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et al. (2002) assume that none of the reference stars have companions and that the measurement 
uncertainties are Gaussian and uncorrelated. Additionally, we assume that multiple planet systems 
do not reduce SIM's efficiency for detecting planets and measuring orbits. Finally, we allow for 
detections of planets with orbital periods more than twice the mission lifetime for which it may be 
difficult to verify that the astrometric deviations from a model with no planets are Keplerian or that 
they are due to the gravitational perturbation of a planetary companion. One change that could 
significantly increase the number of planets SIM would find would be that the planetary mass 
function is much larger for terrestrial-mass planets than implied by extrapolating the planetary 
mass function for giant planets. 



3. Results 

We now present the results of our Monte Carlo simulations. In addition to discussing several 
possible sets of mission parameters, we frequently divide our results into three categories based on 
the planet mass. Planets with m < 3M e would be extremely interesting since they have a mass 
comparable to the Earth and are likely to be rocky bodies formed by a mechanism similar to that 
which formed the Earth. This category would include all the terrestrial planets in our solar system. 
We also consider planets with m < 20M^ that planet formation theories suggest might still be 
mainly rocky bodies that have not undergone substantial growth via gas accretion. Note that this 
category would include Uranus and Neptune as well as the terrestrial planets in our solar system. 
Finally, we also present results for planets of all masses, but these are still limited to m < WMj 
due to our choice of the planet mass function. This category would include all the planets in our 
solar system as well as 93% of the presently known extrasolar planets. 

In table 2 we present the number of planets which are detected or characterized in our sim- 
ulations for all planets (top), planets with mass less than 20M^ (middle), and planets with mass 
less than 3M ffi (bottom). The different columns provide the results of simulations with various 
mission parameters. The ranges provided represent 90% confidence intervals both in this table and 
throughout this paper. These confidence intervals reflect only statistical uncertainties due to the 
random realization of the mass-period distribution (1), not uncertainties in the distribution itself. 

In Fig. 2 we show the likelihood of finding a given number of planets. The different columns 
consider planets with different maximum masses. The top four rows are for a SIM survey with 
1 fias single measurement precision. The top row is for detecting a planet, the second row is for 
measuring the mass and orbital parameters with 30% accuracy, the third row is for measuring the 
masses with 10% accuracy, and the fourth row is for measuring the masses and orbital parameters 
with 10% accuracy. The bottom two rows are for a radial velocity survey of the same stars with 3 
m s^ 1 single measurement precision. The next to bottom row is for detecting the planet, and the 
bottom row is for measuring the masses and orbital parameters. The solid lines are for five-year 
surveys and the dotted lines are for 10 year surveys. We conclude that a SIM survey of 120 stars 
at l^as precision is expected to detect ~ 24 ± 7 planets over five years or ~ 33 ± 8 planets over 
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ten years. The five-year survey would measure the masses and orbits of ~ 16 ± 6 of these planets 
with 30% accuracy and ~ 13 ± 6 planets with 10% accuracy. However, most of these planets are 
relatively massive. If we restrict ourselves to planets with mass < 3M^, then even the ten-year 
survey would only detect ~ 5 ± 3 planets and would probably not measure any masses or orbits to 
10% accuracy (see Table 2). 

In Fig. 3 we explore the effect of varying the mission parameters. The solid line is for a mission 
with 1 /uas single measurement accuracy, the long dashed line for 1.4 fias, the dotted line for 2 ^as. 
The thick lines are for a five-year mission and the thin lines are for a ten-year mission. The upper 
panel shows the number of planets detected and the other panels show the number of planets for 
which masses and orbits are determined. The left panels are for all planets, and the right panels are 
for planets with m < 20M ffi . Note that more planets of all masses are found by missions targeting 
a larger number of stars at lower precision; however, such missions are much less sensitive to the 
low-mass terrestrial planets that may be most interesting. 

In Figs. 4 and 5 we plot the cumulative number of planets found below a given mass. In Fig. 
4 the different line styles are for detection, estimating the mass and orbital parameters with 30% 
accuracy, measuring the mass with 10% accuracy, and measuring the mass and orbital parameters 
with 10% accuracy, while in Fig. 5 the line styles are for surveys with different single measurement 
accuracies and hence different numbers of target stars. The most effective survey for planets less 
than a given mass can be determined by which curve has the largest value for a given mass in 
Figure 5. For example, with a five-year baseline, observing 480 stars with 2 fias precision would 
maximize the total number of planets whose masses are measured with 10% accuracy, but all three 
strategies have a similar yield of planets with mass < 20M e whose masses would be measured with 
10% accuracy. 

In Fig. 6 we plot the cumulative number of planets found below a given mass, divided up by 
the spectral type of the star of each planet. Note that K and M stars are typically closer and lower 
mass than F and G stars, increasing the number of planets found around K and M stars relative to 
F and G stars. This effect is more significant for smaller planetary masses, particularly for masses 
< 10M®. Note that our assumption that the planet mass function is given by equation (1) rather 
than equation (2) results in a larger number of detectable low mass planets around these stars in 
our simulations. 

In Figs. 7 and 8 we plot realizations of our Monte Carlo simulations to illustrate the periods 
and masses that would be found for various mission parameters. In Fig. 7 we require only detecting 
the planet. Of the planets which are detected in our simulations, roughly half have an orbital period 
longer than the mission lifetime and ~ 40% have an orbital period more than twice as long as the 
mission lifetime (for both five- and ten-year missions). The long lever arm of large orbits partly 
compensates for the reduction in sensitivity due to the larger orbital period. In Fig. 8 we only plot 
planets for which masses and orbits are measured with 10% accuracy. These figures illustrate that 
during a five-year mission SIM would be able to detect planets in the mass range 0.01 — O.lMj with 
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orbital periods of several years that are inaccessible to current radial velocity surveys. However, 
SIM is far more effective when it has the longest possible mission lifetime. A ten-year mission is 
far more likely to detect planets that cannot be detected by current radial velocity surveys. Also 
note that there is considerable variation between realizations, particularly for lower mass planets. 

In Fig. 9 and Fig. 10, we show the number of planets that a two-tier survey similar to EPIcS 
would likely discover and the cumulative distribution of the masses such a survey would likely 
detect. During a five- or ten-year mission, our two-tier observing strategy would be expected to 
detect ~ 108 ± 21 or 162 ± 24 planets and measure masses and orbits with 10% accuracy for 
~ 50 ± 15 or 84 ± 18 of these, more planets than would be detected by any of the other SIM 
observing strategies that we examined. During a five or ten year mission, a two-tier survey would 
be capable of detecting a planet with a mass ~ 1.5 or 0.7 M® and measuring the mass and orbital 
parameters with 10% accuracy for a planet with mass ~ 8M§ or ~ Thus, a two-tier observing 

strategy similar to EPIcS could allow SIM to discover and characterize numerous giant planets, 
while maintaining a (somewhat reduced) possibility of detecting terrestrial-mass planets. 

4. Discussion 

In this section we consider several questions about the planet finding capabilities of SIM in 
light of our simulations. 

4.1. Will SIM detect Earth-mass planets? 

Of the mission parameters that we consider, a program with 1 /ias single measurement accuracy 
and 120 target stars is most effective for detecting the lowest mass planets (see Figs. 4 and 5). With 
these parameters, it is likely that SIM would detect one or two planets with m < 1M© in a ten-year 
mission, respectively. While a ten-year mission might measure the mass and orbital parameters of 
one 1M® planet with 30% accuracy, either mission would be extremely unlikely to measure either 
the mass or orbital parameters with 10% accuracy. 

If the criterion for "Earth-mass planets" is relaxed to m < 3M e , then ~ 2 ± 2 or ~ 5I3 
(90% confidence intervals) such planets would be expected to be discovered from a five or ten-year 
mission. Masses and orbital parameters might be measured with 30% accuracy for 1 ± 1 or 2±% 
such planets and with 10% accuracy for +2 or 1^ such planets. The two-tier strategy that we 
consider would be nearly as effective as a 120 star 1 ^uas survey for detecting and measuring orbital 
parameters of planets with masses < 3M© (see Fig. 9). 

It is important to note once again that in this mass range, we are relying on a significant 
extrapolation of the observed planetary mass function. If the actual frequency of Earth-mass 
planets is significantly greater than suggested by present radial velocity observations, then the 
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number of detected Earth-mass planets could be higher. If all stars have an Earth-mass planet at 
1 AU, a ten-year survey of 120 stars at 1 fiss precision could detect ~ 10 ± 4 planets, and measure 
the masses and orbital parameters with 30% accuracy for ~ 3 ± 2 of these. Further, SIM would 
measure the mass and orbital parameters to within 10% for 1+^ such planets. 

4.2. Will SIM detect terrestrial-mass planets? 

Traditional theories of planetary formation (Safronov 1969) predict that the giant planets form 
by accretion of gas onto rocky cores of 5 — 20M e , while lower mass planets form via a different 
mechanism. Here we evaluate the likelihood for SIM to detect planets with m < 20M e . 

If the primary goal is to detect planets with m < 20M e , then of the four observing plans we 
have considered, the one with 2 fias single measurement precision is optimal (see Fig. 5a) and could 
be expected to detect 11 ± 6 or 25 ± 9 planets with masses < 20M® for five- and ten-year missions, 
respectively (see Fig. 3 top). A five- or ten- year mission is expected to make 4+g or 12 ± 6 mass 
and orbital parameter determinations with 30% accuracy (see Fig. 3 middle right) and 3 ± 3 or 8^4 
mass and orbital parameter determinations with 10% accuracy (see Fig. 3 bottom right). 

The two-tier strategy would detect ~ 5^g or ~ 17 ± 9 planets with masses m < 2OM0 during 
a five- or ten-year mission, respectively. It would measure masses and orbits with 30% accuracy 
for 2^2 or planets and measure the masses and orbital parameters with 10% accuracy for 2+_l 
or 5^4 such planets depending on the length of the mission. 

Again, it is important to remember that these estimates depend on an extrapolation of the 
planet mass function observed for larger masses, although in this case the extrapolation in mass is 
relatively small, since radial velocity surveys already detect ~ 30M e planets. 

4.3. What is the smallest mass planet that SIM is likely to detect? 

Any of the ten-year SIM surveys which we considered is likely to detect one planet with a mass 
~ 1M®, for any of the target lists we consider. Similarly, a 1 /ias five-year SIM survey is likely 
to detect one planet with a mass ~ I.5M0 and the other five-year surveys would likely detect one 
planet with a mass ~ 2 — 3M ffi (see Figs. 4 and 10). Of course, the mass of the smallest planet that 
SIM will actually detect may be significantly larger or smaller than these predictions due to small 
number statistics and uncertainties in the mass function. 
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4.4. What is the smallest mass planet that SIM is likely to get orbits for? 

Either a 1, 1.4, or 2 /xas survey is likely to measure with 30% accuracy the mass and orbital 
parameters of one planet with a mass ~ 5M e for a five-year survey. During a ten-year survey, SIM 
would likely measure with 30% accuracy the mass and orbital parameters of one planet with a mass 
~ 1.5M® with the 1 ^as survey and ~ 2M e with either the 1.4 or 2 /xas surveys. The two-tier 
strategy would likely measure with 30% accuracy the mass and orbital parameters of one planet 
with a mass ~ 5 or 1.5Mq for five- or ten-year surveys. 

If 10% accuracy is required, then a 1 /was survey is likely to measure the mass and orbital 
parameters of one planet with a mass ~ 7 or 2M® for five- or ten-year surveys, respectively. 
Similarly, a 2 /ias survey is likely to measure with 10% accuracy the mass and orbital parameters 
of one planet with a mass ~ 8 or 3M e for five- or ten-year surveys, respectively (see Fig. 4). The 
two-tier strategy would likely measure the mass and orbital parameters with 10% accuracy for a 
~ 8M e planet in five years and a ~ 3M e planet in ten years (see Fig. 10). 

4.5. How many new planets will SIM detect? 

We can estimate the number of planets of all masses that SIM is likely to find. These estimates 
are much more reliable than the previous ones, since radial velocity surveys are sensitive to the 
typical planets that SIM is expected to find. Our estimates for a five-year mission duration are also 
more reliable than for a ten-year mission, since radial velocity surveys are just beginning to detect 
planets with orbital periods ~ 10 years. 

If its primary goal is to detect and measure masses and orbits for the greatest number of 
planets, regardless of mass, then SIM should pursue a survey of a large number of stars with a 
relatively low single measurement accuracy. For example, using a single measurement precision of 
1 or 2 fias, SIM would be expected to detect ~ 27 ± 7 or 69 ± 13 planets in a five-year mission or 
~ 33 ±8 or 98 ±15 planets in a ten-year mission (see Fig. 4). The two-tier strategy that we consider 
would likely discover ~ 108 ± 21 or 162 ± 24 planets for a five- or ten-year mission, respectively 
(see Fig. 9). 

Radial velocity surveys are also capable of discovering giant planets with orbital periods of 
several years. Thus, it is important that we estimate the number of "new" planets which SIM will 
detect that would not be discovered by a radial velocity survey of the same stars. A ten-year radial 
velocity survey with 3 m s" 1 precision would be expected to detect ~ 8 ± 4 and 34 ± 9 planets, for 
the 1 and 2 fias target lists. For the same two target lists, SIM would be expected to detect ~ 17 
or 44 new planets in a five-year mission and ~ 25 or 68 new planets in a ten-year mission (see Fig. 
4). The two-tier strategy would be expected to detect ~ 57 or 97 new planets for a five- or ten-year 
mission. However, conducting such a large radial velocity survey and including a significant fraction 
of M stars would require a large amount of observing time. 
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4.6. How many new planets will SIM get orbits for? 

Unfortunately, SIM will not obtain accurate planetary masses or orbits for many of the planets 
it detects. For example, using a single measurement precision of 1 or 2 //as, SIM would be expected 
to determine masses and orbits with 30% accuracy for ~ 16 ± 6 or 44 ± 11 planets for a five-year 
mission or ~ 22 ± 8 or 66 ± 12 planets for a ten-year mission (~ 65% of the planets detected with 
the same mission parameters). SIM would be expected to determine masses and orbits with 10% 
accuracy for ~ 13 ± 6 (54%) or 36 ± 9 (52%) planets for a five-year mission or ~ 19 ± 7 (58%) or 
55 ± 11 (56%) planets for a ten-year mission, depending on the single measurement precision. 

Ten-year radial velocity surveys of the same stars with 3 m s _1 accuracy would determine 
accurate orbits (but not m and sini separately) for ~ 6 ± 4 or 25 ± 8 planets, for the 1 and 2 /xas 
target lists. Depending on the single measurement accuracy, SIM would be expected to determine 
masses and orbits with 30% accuracy for ~ 10 or 26 new planets (i.e., planets without orbits 
determined by the radial velocity survey) for a five year mission and ~ 16 or 44 new planets for 
a ten-year mission. The two-tier strategy could measure masses and orbits with 30% accuracy for 
~ 66 ± 16 and 106 ± 20 planets for a five or ten year mission, respectively, of which ~ 31 and 60 
would not be measured by a radial velocity survey of the same stars. Of course, radial velocity 
measurements alone cannot determine the mass and inclination separately. 

4.7. What is the effect of increasing the number of target stars? 

If additional SIM observing time were to be allocated to planet surveys, then SIM would be 
expected to detect a larger number of planets. However, the additional target stars would tend to 
be more distant than the stars already targeted, so the increase is less than linear in the number 
of target stars. In Fig. 11 we plot the number of planets detected (solid lines) as a function of 
the number of target stars, assuming an equal number of F, G, K, and M stars are targeted and 
ad = Ifias. The dash-dotted lines are for measurements of the mass and orbital parameters with 
30% accuracy, the dashed lines are for measurements of the mass with 10% accuracy, and the 
dotted lines are for measurements of the mass and orbital parameters with 10% accuracy. The top 
row shows the number of detections of planets of all masses, the middle row shows the number of 
detections of planets less massive than 20M®, and the bottom row shows the number of detections 
of planets with masses less than The left column is for a five year mission with 24 two 

dimensional observations of each target star, and the right column is for a ten year mission with 
48 two dimensional observations of each target star. Clearly, searches for Earth- mass planets could 
greatly benefit from a significant increase in the number of stars targeted. In principle, a ten 
year mission could be subdivided into two five year surveys with separate target lists. While two 
separate five year surveys would detect a larger total number of planets (by ~ 35%), a single ten 
year survey would be expected to detect slightly more planets with masses less than 20M e (by 
~ 10%) and a few times as many planets with masses less than 3M e . Therefore, if the SIM mission 
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can be extended beyond the planned five years, planet searches should continue observing stars 
which were part of the original target list. 

An alternative strategy for a ten-year mission is to double the number of target stars and 
still make only 24 two dimensional observations of each target star. This would result in an even 
greater number of planets detected and characterized to within 30% for both planets of all masses 
(by ~ 70%) and planets with mass less than 20M® (by ~ 35%). The average number of planets 
characterized with masses less than 3M e would be reduced by ~ 10% compared to the a ten year 
survey with twice as many observations of half as many target stars. 

4.8. How many radial-velocity planets will SIM get orbits for? 

If SIM were to observe each extrasolar planet already known from radial velocity surveys, 
then it would be able to make 10% accurate measurements of the mass and orbital elements for 
80 ±3, 74 ± 4, or 67 ± 5 of the 99 known planets using 24 two dimensional measurements at 1, 
2, or 4 juas precision (see Fig. 12). Of the eleven known multiple planet systems, 1 //as or 2 //as 
astrometry would allow accurate orbital determinations for at least two planets in seven or five of 
those systems, assuming that SIM's ability to measure orbits is not diminished in multiple planet 
systems. In practice, the primary importance of these measurements would be to confirm the 
interpretation of the radial velocity measurements and to constrain the inclinations of the orbits. 

5. Conclusions 

All of our results are based on the assumption that the mass function for terrestrial-mass 
planets can be obtained by extrapolating the mass function for giant planets determined by radial 
velocity surveys. If this assumption is correct, then we find that a ten-year planet search with SIM 
might detect a few Earth mass planets. SIM could be expected to measure the masses and orbital 
parameters with 10% accuracy for planets with masses as small as ~ 2M e . While SIM could find 
more low mass planets if they are more common than current estimates, a null result would not 
demonstrate that planets with masses less than ~ 10M^ are less common than presently expected 
on the basis of extrapolating the results of radial velocity surveys. Including giant planets, SIM is 
expected to detect ~ 24 — 162 planets; the higher number results from the ten-year surveys which 
target a large number of stars at relatively low precision. Of these, a significant fraction, 25 — 65%, 
could also be discovered by a 3 m s^ 1 radial velocity survey. With the observing strategies that we 
have examined SIM will measure the masses and orbits of ~ 16 — 106 planets with 30% accuracy 
and the masses and orbits of ~ 13 — 84 planets with 10% accuracy. Of these, ~ 30 — 70% could be 
measured by a 3 m s _1 radial velocity survey, 
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Distance (pc) Magnitude (V) 



Fig. 1. — Cumulative number of stars in our target lists as a function of distance (left) and 
magnitude (right). The top, middle, and bottom rows are for the target lists with 120 stars (ad = 1 
fias), 240 stars (o^ =1.4 fias), and 480 stars (cr^ = 2 /xas). The long dashed, dotted, short dashed, 
and dashed-dotted lines are for F, G, K, and M stars, respectively. The solid bold lines are for all 
the stars in each list (the normalization of this line has been divided by four). 
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Fig. 2. — The likelihood of finding a given number of planets. The top four rows are for a SIM 
survey with a single measurement accuracy of 1 /xas that targets 120 F, G, K, and M stars (30 in 
each spectral class). The top row is for detections, the second row is for estimates of the mass and 
orbital parameters with 30% accuracy, the third row is for mass measurements with 10% accuracy, 
and the fourth row is for measurements of the mass and orbital parameters with 10% accuracy. The 
bottom two rows are for a radial velocity survey of the same stars with 3 m s _1 single measurement 
accuracy. The fifth row is detections and the bottom row is orbital determinations. The solid lines 
are for five-year surveys and the dotted lines are for ten- year surveys. In some panels the two lines 
coincide. The left column is for planets of all masses, the center column is for planets with mass 
less than 20M e , and the right column is for planets with mass less than 3M®. 
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Fig. 3. — The likelihood of finding a given number of planets. In each panel, the solid line is for a 
mission with 1 fias single measurement accuracy, the long dashed line for 1.4 fias, the dotted line 
for 2 //as. The thick lines are for a five-year mission and the thin lines are for a ten-year mission. 
The top panels are for detecting planets, the middle panels are for measuring the mass and orbital 
parameters with 30% accuracy, and the bottom panels are for measuring the mass and orbital 
parameters with 10% accuracy. The left panels are for all planets, while the right panels are for 
planets with mass less than 20M e . 
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Fig. 4. — The results of Monte Carlo simulations to estimate how many planets would be found 
below a given mass. We present the weighted average of thousands of realizations to yield accurate 
statistics. The solid lines are for detections, the dash-dotted lines are for measurements of the mass 
and orbital parameters with 30% accuracy, the dashed lines are for measurements of the mass with 
10% accuracy, and the dotted lines are for measurements of the mass and orbital parameters with 
10% accuracy. Each panel is for a different set of mission parameters and target lists. The top 
panels (a and b) are for 1 fiss single measurement accuracy and 120 target stars. The middle panels 
(c and d) are for 1.4 fias single measurement accuracy and 240 target stars. Bottom panels (e and 
f) are for 2 /ias single measurement accuracy and 480 target stars. The left panels (a, c, and e) are 
for a five-year mission, while the right panels (b, d, and f) are for a ten-year mission. The masses 
of Jupiter and the Earth are indicated by J and ©, respectively. 



-20- 



100 



10 



1 



i i iiiiii i i i inn 


inn i i 1 1 mi 


— a 




- Detections 




- / / 






J 


w/>\ 


t 



10 



1 



TTTT 
b 



TTTI 



ttti 



I I I 1 1 111! 




10 _ 6 1Q _5 10 _4 1Q -3 1Q -2 1Q-6 1Q -5 1Q -4 1Q -3 1Q -2 



I 100 



10 



I I I II 1 I I I INI 

c 

Masses (10%) 



I I I 14 




10 



1 

10 _ 5 1Q -4 1Q -3 10-2 1Q-5 

Mass (M ) 



1 IIIIII 1 1 1 1 1 III 1 

: d 


1 1 1 IILL 


" Orbits (10%)...--" 








-/// 




if j 




/■ 1 1 






io- 



10 



10 



-2 



Fig. 5. — The results of Monte Carlo simulations to estimate how many planets would be found 
below a given mass. We present the weighted average of thousands of realizations to yield accurate 
statistics. The solid, long dashed, and dotted lines are for missions with a single measurement 
accuracy of 1 ;uas for 120 target stars, 1.4 //as for 240 stars, and 2 //as for 480 stars, respectively. 
The thick lines are for a five-year mission and the thin lines for a ten-year mission. Panel a is for 
a detection, panel b for a 30% measurement of the mass and orbital parameters, panel c for a 10% 
measurement of the mass, and panel d for a 10% measurement of the mass and orbital parameters. 
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Fig. 6. — The results of Monte Carlo simulations to estimate how many planets would be found 
below a given mass. We present the weighted average of thousands of realizations to yield accurate 
statistics broken down according to the spectral type of the star. The long dashed lines are for 
F stars, the dotted lines are for G stars, the short dashed lines are for K stars, and the dashed- 
dotted lines are for M stars. Each panel is for a different set of mission parameters and target 
lists. The top panels (a and b) are for a 1 fiss single measurement accuracy and 120 target stars 
equally distributed among the four spectral types. The middle panels (c and d) are for a 1.4 //as 
single measurement accuracy and 240 target stars. Bottom panels (e and f) are for a 2 /tas single 
measurement accuracy and 480 target stars. The left panels (a, c, and e) are for a five-year mission, 
while the right panels (b, d, and f) are for a ten-year mission. 
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Fig. 7. — In each panel we show the masses and periods of planets detected from a single simulation. 
The left and right columns are for five- and ten-year missions, respectively. The first, second, and 
third rows are for 1 //as, 1.4 //as, and 2 //as single measurement accuracy, respectively. The bottom 
row is for a two-tier target list, including 52 stars at 1 //as precision and 1092 stars at 4 //as 
precision. In each panel, open squares and solid squares are for planets detected by SIM only and 
by both SIM and radial velocity searches. The © and J symbols indicate the masses and orbital 
periods of the Earth and Jupiter. 
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Fig. 8. — In this figure we show the masses and periods of planets whose masses and orbital elements 
are measured with 10% accuracy. A single simulation is shown in each panel. The left and right 
columns are for five- and ten-year missions, respectively. The first, second, and third rows are for 1 
yuas, 1.4 /xas, and 2 //as single measurement accuracy, respectively. The bottom row is for a two-tier 
target list, including 52 stars at 1 //as precision and 1092 stars at 4 //as precision. In each panel, 
open squares and solid squares are for planets detected by SIM only and by both SIM and radial 
velocity searches. The © and J symbols indicate the masses and orbital periods of the Earth and 
Jupiter. 
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Fig. 9. — Here we show the likelihood of finding a given number of planets using a two-tier strategy. 
The top four rows are for a SIM survey that targets the 52 nearest stars with a single measurement 
accuracy of 1 ,uas and 1092 conveniently located nearby stars with a single measurement accuracy 
of 4 /xas. The top row is for detections, the second row is for estimates of masses and orbital 
parameters with 30% accuracy, the third row is for mass measurements with 10% accuracy, and 
the fourth row is for measurements of the masses and orbital parameters with 10% accuracy. The 
bottom two rows are for a radial velocity survey of the same stars with 3 m s _1 single measurement 
accuracy. The fifth row is detections and the bottom row is orbital determinations. The solid lines 
are for five-year surveys and the dotted lines are for 10 year surveys. The left column is for planets 
of all masses, the center column is for planets with mass less than 20M®, and the right column is 
for planets with mass less than 3M®. 
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Fig. 10. — Here we show the results of a Monte Carlo simulation to estimate how many planets 
would be found as a function of maximum mass using a two-tier strategy. We present the weighted 
average of thousands of realizations to yield accurate statistics. The solid lines are for detections, 
the dash-dotted lines are for 30% estimates of orbital parameters, the dashed lines are for 10% mass 
determinations, and the dotted lines are for 10% orbital determinations, The thick lines are for a 
five-year mission and the thin lines for a ten-year mission. These results are for a SIM survey that 
targets the 52 nearest stars with a single measurement accuracy of 1 fias and 1092 conveniently 
located nearby stars with a single measurement accuracy of 4 //as. The masses of Jupiter and the 
Earth are indicated by a J and ©, respectively. 
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Fig. 11. — In this figure we plot the number of planets detected as a function of the number of 
target stars, assuming an equal number of F, G, K, and M stars are targeted and ad = 1/xas. 
The solid lines are for detections, the dash-dotted lines are for measurements of the mass and 
orbital parameters with 30% accuracy, the dashed lines are for measurements of the mass with 10% 
accuracy, and the dotted lines are for measurements of the orbital parameters with 10% accuracy. 
The top row shows the number of detections of planets of all masses, the middle row shows the 
number of detections of planets less massive than 20M®, and the bottom row shows the number 
of detections of planets with masses less than 3.M®. The left column of panels is for a five year 
mission with 24 two dimensional observations, and the right column of panels is for a ten year 
mission with 48 two dimensional observations. Note that increasing the number of target stars 
above ~ 300 in a five-year survey does not substantially increase the number of < 3M ffi planets 
detected or characterized. 
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Fig. 12. — In this figure we show the number of known extrasolar planets with a given astrometric 
signal. The vertical lines show the minimum astrometric signature necessary for SIM to measure 
the mass and orbital elements to within 10% with 95% probability according to the simulations 
of Sozzetti et al. (2002). The dashed, dotted, and dashed-dotted lines assume single measurement 
precisions of 1, 2, and 4 //as, respectively. We have assumed random inclinations of the orbital 
plane to the sky and averaged several realizations, resulting in bins that have non-integer heights. 
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Table 1. Parameters for detection and measurement probability functions 

Orbital Period a Detection 13 30% Mass & Orbit c 10% Mass c 10% Orbit c 

0—0.2 0.77 1.77 2.81 3.11 

1.82 2.77 3.21 5.01 

2.78 1.95 8.88 9.15 

5.59 7.09 1.95 2.14 

0.2—0.8 0.86 2.03 3.10 1.24 

1.86 2.97 3.74 4.55 

2.98 1.39 7.09 3.25 

5.34 6.58 1.43 1.24 

0. 8—1 1.73 2.78 3.69 7.21 

1.81 3.48 4.45 8.09 

3.62 1.31 5.69 3.11 

7.11 5.55 1.24 1.30 

1. — 1.4 2.25 4.15 4.34 6.70 

2.29 5.91 8.09 21.9 

5.90 6.77 6.84 4.62 

3.25 1.40 1.01 1.41 

1.4—2 3.08 5.89 6.35 7.04 

3.92 11.5 16.2 47.5 

6.12 6.85 6.31 9.57 
2.40 1.41 0.99 1.30 

2— 3 4.37 9.08 11.5 13.4 

8.31 21.1 26.9 76.2 

7.09 6.06 4.34 5.70 
2.12 1.59 1.22 1.27 

3— 4 11.8 11.7 27.3 42.1 

12.7 37.3 36.0 94.6 

2.10 5.98 2.87 2.87 
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Table 1 — Continued 


Orbital Period a 


Detection 13 


30% Mass & Orbit c 10% Mass c 


10% Orbit c 




2.14 


1.94 1.40 


1.10 


4—8 


23.9 








38.6 








1.75 








1.77 






8—12 


10.8 








145. 








1.88 








1.93 







Note. — Each entry lists the parameters b, c, j3, and 7 (from top to bottom) for use in equation 
(5). 

a Divided by the time span of observations 

b For a 0.1% false alarm rate 
c For 95% of the planets detected 
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Table 2. The number of planets detected or characterized in our simulations 





1 fias 


5 Year 
2 fias 


2 Tier 


1 /xas 


10 Year 
2 /Ltas 


2 Tier 


All Masses 














Detect 


24 ±7 


69 ±13 


108 ± 21 


33 ±8 


98 ± 15 


162 ± 24 


Mass & Orbit (30%) 


16 ±6 


44 ±11 


66 ± 16 


22 ±8 


66 ± 12 


106 ± 20 


Mass & Orbit (10%) 


13 ±6 


36 ±9 


50 ± 15 


19 ±7 


55 ±11 


84 ± 18 


m < 20M e 














Detect 


7±5 


11 ±6 


5t§ 


13 ±6 


25 ±9 


17 ±9 


Mass & Orbit (30%) 


4±3 




2^ 


8±4 


12 ±6 


7 +6 


Mass & Orbit (10%) 


9+3 
z -2 


3±3 


9+3 
z -2 


c;+4 
°-3 


8tl 


r+5 
°-4 


m < 3M e 














Detect 


2±2 






rc+4 
°-3 






Mass & Orbit (30%) 


1 ± 1 


0+ 1 


+2 


9+3 
z -2 




2±2 


Mass & Orbit (10%) 


+2 


0+ 1 


0+2 


11? 


11? 





Note. - 



Error bars correspond to 90% confidence intervals. 



